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Abstract 

The quantum dynamics of a driven single-band tight-binding model 
with different boundary conditions is considered. The relation between the 
Hamiltonian describing the single-band tight-binding dynamics and the 
Hamiltonian of a discrete-charge mesoscopic quantum circuit is elucidated. 
It is shown that the former Hamiltonian, with Dirichlet boundary condi- 
tions, can be considered as a realization of the deformed parafermionic 
polynomial algebras. 
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1 Introduction 

The dynamics of a quantum particle in a periodic potential under the effect of 
an external field is one of the most fascinating phenomena of quantum physics 
[1,2]. Under this condition, the electronic wave function displays so called Bloch 
Oscillations (BO) [3,4], the amplitude of which is proportional to the band 
width. A suitable orthogonal basis for investigating structures with periodic 
potential is using the localized states. These states are also called the Wannier- 
Stark states [5-9] and the nature of this states has a significant influence on the 
electronic transport properties of solids. The electronic BO was observed for the 
first time in semiconductor superlattices [10]. In [10,11], the first experimental 
observation of BO was reported. Recently, an increasing interest in the dynamics 
of BO can be observed [7,11-20]. 
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An explicit time dependence of physical properties can appear under several 
conditions. The most obvious case is when the external field depends on time. 
In principle, one has to define an initial state \^>(t = 0)) and solve the time- 
dependent Schrodinger equation. This is a partial differential equation with at 
least two variables which is separable in some simple cases [21,22], but rarely 
solved analytically. 

In particular for the tight-binding model, some analytical expressions have 
been derived both for time-dependent and -independent fields [23-25] . Recently 
in the context of the tight-binding model, a treatment based on the dynamical 
Lie algebra is proposed in [26-29]. The time-dependent Hamiltonian which we 
will mainly consider here is the Hamiltonian of a finite system with Dirichlct 
boundary conditions defined by 

h = + 1| + \j + i)od + Fwjhjmi (i) 
.?=i .?=i 

where the ket \j) represents a Wannier state located on the site j. The 
Hamiltonian (1) with infinite boundary conditions has been investigated for 
example in [15,27]. The Hamiltonian (1) can be treated as the Hamiltonian of 
a finite quantum wire, in tight-binding approximation, under the influence of 
a time-dependent electric field. This Hamiltonian appears in the discussion of 
the quantum transport through nanoscale conductors [30]. Hamiltonian (1) also 
appears in models describing the dynamics of cold atoms in optical lattices in 
tight-binding approximation [31]. 

In the present work we elucidate the similarity between the Hamiltonian 
(1) with infinite boundary conditions, i.e, — oo < j < oo, and the Hamiltonian 
describing the discrete-charge mesoscopic quantum circuits [32,33,34]. In par- 
ticular, we show that the Hamiltonian (1) with Dirichlct boundary conditions 
can be considered as a realization of the deformed parafcrmionic polynomial 
algebras. Parafermions of order p (with p being a positive integer) have been 
introduced in [35,36,37]. The nature of these particles is such that at most p 
identical particles can be found in the same state. The usual spin half fcrmions 
correspond to p = 1. The notion of parafcrmionic algebra has been extended by 
Qucsnc [38]. The relation between parafermionic algebras and other algebras 
has been given in [39,40,41]. 

Some extensions of the Hamiltonian (1) are in order. For example one can 
add an aperiodic potential to the Hamiltonian , see for example [42] or add non 
liner terms as in pendulum model [43]. 

The layout of the paper is as follows: In section 2, we introduce the model 
and discuss in summary the discrete-charge mesoscopic quantum circuits. In 
section 3, periodic boundary condition is discussed. In section 4, a realization 
of the deformed parafermionic polynomial algebra is presented. In section 5, the 
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time-evolution operator for a finite system with Diriehlet and periodic boundary 
conditions is discussed. In section 6, we present a summary of our results. 



2 The model 



The system under study is the derived quantum motion of a charged particle 
in a one-dimensional array of single-band quantum wells in tight-binding ap- 
proximation and under the action of an arbitrary time-dependent external field 



where the ket \j) represents a Wannier state located on site j. These states 
fulfill the orthogonality condition {j\j'} — 5jj> . The force F(t) is an arbitrary 
time-dependent force and the real parameter G is the nearest-neighbors cou- 
pling strength. In the following we will investigate the Hamiltonian (2) and its 
applications in derived discrete-charge quantum circuits and specially we will 
show that this Hamiltonian, with closed boundary conditions, can be considered 
as a realization of the generalized parafermionic polynomial algebras. 

2.1 Infinite chain 

As the first boundary condition let us assume that the Hamiltonian (2) describes 
an infinite chain of single-state quantum wells, i.e, — oo < j < +oo. In this case 
the Hamiltonian (2) can be written as [27] 



F(t). 



The Hamiltonian of such a system can be written as 



H = Gj2(\j)(j + 1| + \J + 1)01) + F(t) £ j|j)0'l 



(2) 




H = G(K + K^) + F(t)N, 



(3) 



where 



oo 



K 



j=-oo 



oo 



j=-oo 



oc 



N 




(4) 



Operators K and act as ladder operators on the Wannier states 



K\3) = 1.7-1) 
K ] \j) = |.7 + 1) 



(5) 
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and the operator N acts on the Wannier states as the position operator 



N\j)=j\j). (6) 
These operators fulfill the following algebra 

[N,K] = -K, 

[n,k*] = k\ 

[k\k] = o. (7) 

A quantum theory for mesoscopic circuits was proposed by Li and Chen 
[32] where charge discreteness was considered explicitly. In this case the charge 
operator q takes on discrete eigenvalues when acting on the charge eigenvectors 

q\n)=nq e \n), (8) 

where n € Z (set of integers) and q e is the electron charge. The charge 
operator q can be represented in terms of the localized states as 

oo 

Q = <leJ2 3\j)U\=QeN. (9) 
j=—oo 

Discrete-charge mesoscopic quantum circuits are described by the Hamilto- 
nian 

H = -^(Q + &-2) + V(q), 
i 

Q = eh , 
[q,p] = ih, 

(10) 

where L is the inductance of the circuit. The operators Q and are charge 
ladder operators which together with the charge operator q fulfill the following 
commutation relations 



[q, Q] = -QeQ, 

[q,&)=q e Q\ 

[Q + ,Q] = 0. (11) 
Therefore the operators N = q^ 1 !}, Q and satisfy the same algebra (7). 

-2 

When V(q) — that is when there is a capacitor in the circuit with capacity 
C, we have a LC-design and for V(q) = 0, we have a pure L-design. 
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The timc-dependent Hamiltonian of a L-design circuit in charge representa- 
tion and under the influence of an external potential e(t) is [32] 

H=-^(Q + &) + e(t)q, (12) 
which is mathematically equivalent to the Hamiltonian (3). 

3 Periodic boundary condition 

We can also take the periodic boundary condition i.c, \n + N) = \n). This kind 
of boundary condition appears for example in describing a circular quantum 
wire in tight-binding approximation. Where, a time-dependent magnetic field 
piercing the ring, acts as a driving force according to the Faraday's law. In this 
case the Hamiltonian (2) and the corresponding operators are defined as 

N N 

H = GJ2(\J)U + 1| + \J + 1)01) + F(t) £ ili)01, (13) 
i=i i=i 



N 



& = Eli + i)0'l, 



N 

n = x>->oi- (14) 

The commutation relations are now a deformation of the algebra (7) as 



[N,K] = -K{1 - h(N)), 
[N,K*] = (l-h(N))K\ 

[k,k^] = o, 

(15) 



where the polynomial operator h(N) is the deformation operator defined by 

at(--\\n-i i= N 
= (N-l)\ II & - 3) = ^Il)(l|- (I 6 ) 
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4 Deformed parafermionic algebra 



In this section we show that the Hamiltonian (1), with Dirichlct boundary con- 
ditions, can be considered as a realization of the deformed parafermionic poly- 
nomial algebra. The Hamiltonian (1) appears for example in the dynamics of 
cold atoms in optical lattices [31,44] and also in the modeling of a finite quan- 
tum wire, in tight-binding approximation and under an external time-dependent 
electric field [30]. The Hamiltonian can be written as 

JV-l N 

h = gJ2 (li)O' + il + li + i)0'l) + n*)E^')oi, 

N-1 
N-1 

& = + i)01, 

jV = f>-)0l- (17) 

The Hamiltonian in terms of the generators can be rewritten as 

H(t) = G(K + Ki) + F(t)N, (18) 
the new commutation relations arc now as follows 



[N,K] = -K, 

[K^,K] = \N)(N\-\1)(1\. (19) 

From (17), it is clear that the operators K and are nilpotent of order N, 
that is, 

K N = K ]N = 0. (20) 

The last commutation relation in (19) can be expressed as a polynomial in 
terms of the operator N. This polynomial is not necessarily unique and we can 
choose a polynomial with the smallest degree. For this purpose let the unknown 
polynomial be f(N) then 

f(N) = \N)(N\-\l)(l\. (21) 



G 



The polynomial f(N) can be determined from (21) using the following con- 
ditions 



f(N)\n) = 0, for 2 < n < N - 1, 

f(N)\N) = \N), 
f(N)\l) = 

after some simple algebra we find 



(22) 



N-l 

W) = Jj^q[Il("-j)> foroddiV>3, 
1 

/(AO = (N _ lV VN-(N + l))l[(N-j), for even N > 2, 
f(N) = (27V -3) foriV = 2. 

(23) 

Let us define the polynomial operator g(N) as 

f(N)=g(N+l)-g(N). (24) 
From (23) one can find g(N) up to a constant as follows 

1 N 

9(N) = j[(N - j), for odd N > 3, 

2 W + 2 W 
fl(^) = (^--j^-)II^ -■?')' for even N >2. 

(25) 

Putting everything together, we arrive at the following polynomial algebra 
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[n,k*] = k\ 
[n, k] = -k, 

[k\k]=f(N)=g(N + l)-g(N), 

[k^k,kk^\ = o, 
k N = k^ N = o, 

K^k = l-g(N) =: <f>(N), 

KK^ = 1 - g(N + 1) =: <j>(N + 1), 

(26) 

where 

0(1)=^ + 1) = O, (27) 

as is clear from the definition of the operators k and . Now we show that 
the polynomial algebra (26) is in fact a deformation of the usual parafcrmionic 
polynomial algebra C — {B, £?t , M} [35] defined by 

[M, B] = -B, 
[M,&] = B\ 
B p+1 = S tp+1 = 0, 

= M(p+1 — M) = [M], 
B& = (M + l)(p - M) = [M + l], 

M=\{[B\B]+p), 

(28) 

where p is a nonzero positive integer. For this purpose let us define the 
operators B, St and M as 

B = tJn(N-N)K, 

St = K1\Jn(N - N), 
M = N-l, 

(29) 

then it can be easily shown that the operators B, fit and M, fulfill the 
same algebra (28). It is interesting that for the special cases N — 2,3 the 
algebra (26) coincides with (28) and for N > 4, it is a deformmation of the 
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algebra (28) with the deformation factor y N(N — N) [38]. The deformation 

factor is not invcrtiblc since \J N(N — N)\N) = but we can make it invert- 
ible by introducing an infinitesimal parameter e and defining a new operator 

D e := \J N(N + e — N) = \J (M + 1)(N — 1 + e — M) which tends to the de- 
formation factor when e — > 0. Now the Hamiltonian can be written in terms of 
the generators of the parafermionic polynomial algebra (28) as 

H = G{t)- l {M)B + Bt^-i(M)) + F(t)(M + 1). (30) 

The propagator for the non deformed Hamiltonian 

H* = G(B + B*<) + F(t)(M + l), (31) 

can be obtained exactly from the dynamical lie algebra method [45]. 

4.1 The su(2) representation of the special cases iV = 2,3 

Setting N = 2 in equations (26), ( 23) we arrive at 



[N,K] = -K, 
[N,fr]=K*, 
[K\k] =27V-3, 

k 2 = kv = o, 
n = i + k^k. 

(32) 

From {K, K^} = I and K 2 = K^ 2 = 0, it is clear that the set {K, k\ 1} is 
a realization of the usual spin-half fermionic algebra. By defining 



Jo 

J+ 
J- 



*4 

-k. 



(33) 



we find 



[Jo, J±] = ±4, [4, J-] - 2J , (34) 

which is the usual su(2) algebra. In terms of J± and Jo, the Hamiltonian 
(18) can be rewritten as 
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H(t) = G(J + + J_) + F(t)(J -^). (35) 

For this Hamiltonian the Casimir operator J 2 = J + |(J+J_ + J-J+) is 
a constant of motion. 

Setting N — 3 in equations (26), ( 23) we find 



[N,k] = -k, 

[N, K^] = K\ 
[K\k] =N-2, 

k 3 = k^ = o, 

n = i + k^k + k^k 2 , 

(36) 



again we can define the su(2) generators as 



J+ 


= V2K\ 


k 


= V2K, 


Jo 


= N-2. 



(37) 

The Hamiltonian in terms of this generators becomes 

H{t) = -j=(J + + J-)+F(t)(j -2). (38) 

For TV = 4, 5, f(N) is a cubic polynomial operator (like Higgs algebra) and in 
general, for TV = n, n + 1, f(N) is a polynomial operator of degree n — 1. 



5 Time- evolution operator 

For N = 2, 3 Hamiltonian has su(2) symmetry and the time-evolution operator 
can be obtained exactly. In this case the time-evolution operator U(t) can be 
written as 

U(t) = e-icM^i/W^-lsWj+e^^j-. (39) 
Inserting the above expression in the following evolution equation 

ihjU(t)=HU{t), (40) 
we can find U(t) straightforwardly [45]. 
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For N > 4, the previous method is not successful and we find a perturbative 
expansion for U (t) in this case. For this purpose let us write the time-evolution 
operator as 

U{t) = e -* v ^V(t), (41) 
where ip(t) = f F(t')dt'. Substituting Eq.(41) into the evolution equation 

ihjU(t) = HU(t), (42) 
we find the following equation for V(t) 

ihjV(t) = G{e-rMt)k + eM*W)V(t), (43) 
where the following relations have been used 

(44) 

The relations (44) can be obtained directly from the Bakcr-Hausdorf formula. 
The Eq.(43) can be written in terms of the state vector \ip(t)) = V(t)\ip(0)) as 

= G(e~^K + eM*W)\il>(t)). (45) 

The new Hamiltonian 

G{e-W>K + eM*W), (46) 
has the following eigenvalues and eigenvectors respectively 

A m = 2Gcoso; m , m = 1,2- • • N, 

N 

\\ m , t)=d m J2 e mip(t) Bin(nw m )|n), (47) 

71=1 

where co m = "" r 1 and d m = . 1 , is a normalization coefficient. 

+ VJ2"=i sin 2 n(w m ) 

Therefore, the eigenvalues are independent from the driving force F(t) but the 
eigenvectors are time dependent and make an orthonormal complete basis at 
any time. Now let us expand the state vector as 

JV 

i=i 

Inserting (48) into (45) we find 

N 

ihC e (t) =fvpYj A u' C i'^) + Mt)Ct(t). (49) 
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where 



Z)n=i 71 sin(nw£) sin(nw^ ) 



A u , = (X e (t)\-\X/(t)} 



N 



N 



(50) 



a ^ sin 2 (ntj£) ^ sin 2 (no>£') 

\ n— 1 n— 1 



are the elements of a symmetric matrix A. The Eq.(49) can be written in matrix 
form as 

ihd(t) = htpAc(t) + xc(t). (51) 

where 





- c\(t) - 




cos(ct>i) 




















cos(w 2 ) 








c = 




X = 2G 








cos(w 3 ) • 







_ C N (t) _ 













cos(o>at) 



(52) 

By rescaling the parameters we can assume 2G < 1 and since cos(w m ) < 
1, (1 < m < N), so the matrix A has the property lim^oo X k — > 0. Now the 
equation (51) is similar to the Schrodinger equation with the perturbation A. 
Therefore the solution of (51) can be written as a series expansion in terms of 
the matrix A as 



C{t) = C°(t) + XC\t) + X 2 C 2 (t) + X a C 3 (t) + 



(53) 



substituting (53) into (51) and equating the terms of the same order in A we 
find 



ihd°{t) 



hip AC (0), 

c°(t) + fup\- 1 A\c 1 (t), 



The first equation leads to 

C<°>(t) = e -*^ C (o)( ). 
For n > 1, let us define A„ := X~ n AX n , then 

(l + ii pA n )c n (t) = - l -c^- 1 Ht). 



(54) 



(55) 



(56) 
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The equation(56) can be solved using the Green function method and treating 
C^-^it) as the source term for C {n \t) as 

C n (t) = cxpH<^)A„)C™(0) - - / exp(-t(f(t - t')A n )C ( - n - 1) (t')dt' . (57) 

Jo 

As the boundary condition we can assume C(0) = C°(0) and C^(0) = for 
k > 1. So the equation (57) can be simplified as 

C n (t) = - % - [ exp(-iip(t - i')A„)C (n - 1) (t')^', (n > 1). (58) 

"> Jo 

Now having (55), we can find C^ n \t) up to any order recursively. 

6 Conclusions 

The similarity between the Hamiltonian of a single-band tight-binding model 
and the discrete-charge mesoscopic quantum circuits is explicitly shown. It is 
shown that the former Hamiltonian, with Dirichlet boundary conditions, can be 
considered as a realization of the deformed parafermionic polynomial algebras. 
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